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Abstract. We prove that for the Navier Stokes equation with dissipation 
(—A)", where 1 < a < |, and smooth initial data, the HausdorfT dimension 
of the singular set at time of first blow up is at most 5 — 4a. This unifies two 
directions from which one might approach the Clay prize problem. 



1. Introduction 



The Navier Stokes equation with dissipation (—A)" is given by 

— +u-Vu + Vp= -(-A)"-u, 



(1) 



where u is a time-dependent divergence free vector field in R^. One sets the initial 
condition 



where uo{x) G C;?°(R3). 

It is easy to see that classical solutions to this equation on a time interval [0, T] 
satisfy conservation of energy, namely that 



Jo 

The second term on the right is called the dissipation term. 

Caffarelli, Kohn, and Nirenberg 1^ showed that when a — 1, the singular set of 
a generalized weak solution (a weak solution satisfying a certain generalization of 
conservation of energy) to the system has parabolic Hausdorff dimension 

at most 1. This could be considered a first step towards showing global strong 
solvability. Indeed any improvement in this upper bound on the dimension would 
be genuine progress towards solution of the global solvability problem. 

Another well known fact is the following. (We learned the proof which appears in 
the following section from Diego Cordoba. It appears to exist in the folklore and 
we have not succeeded in attributing it. Sinai and Mattingly |Q recently gave a 
different proof.) 
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u{x, 0) = Uo{x) 



(2) 
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Proposition 1.1. If a > |, one has global strong solvability for the system (^,(||). 



Indeed proposition |L1] could be viewed as a first step towards the solution of the 
global strong solvability for a = 1 and any improvement in the exponent | could 
be viewed as genuine progress. 

The purpose of this paper is to interpolate the two results so that these two paths 
to progress are unified. We prove 

Theorem 1.2. If T is the time of first breakdown for the system (|^),(^), with 
J > a > 1 then the Hausdorff dimension of the singular set at time T is at most 
5 - 4q;. 



Heuristically, the theorem ought to be a mild generalization of Proposition 
We think about the microlocal analysis of a solution u in terms of coefficients uq 
associated to cubes Q. Very roughly speaking, the coefficient uq should be viewed 



as a generalized wavelet coefficient. The proof of Proposition 1.1 goes wrong for 



a < I only because of a small number of cubes with large coefficients. At any 
time, the set of points contained in arbitrarily small such cubes has dimension at 
most 5 — 4a. The difficulty lies in preventing the energy from flowing between far 
away cubes near the time of blow up. This is not quite so easy as it sounds because 
energy can flow from small cubes to large ones containing them and this is a quite 
non-local effect. 



We were also unable to directly generalize the proof of Caffarelli, Kohn, and Niren- 
berg. They rely on what they call the "generalized energy inequality" which is 
built on an amusing property of the divergence free heat equation. Let cj) be any 
compactly bump function in space and time and u a divergence free vector field, 
then 

j ((^ - A)^., ^u)dt = j (^^((«, 0")) + (V^, 0Vu) ~{{\^ + A0)zi, u))dt. 

The first term represents a change in local energy. The second term represents 
local dissipation. The third term is an error which can be made insignificant by 
choosing (p to satisfy the backwards heat equation. It is this method of controlling 
the error which we are unable to generalize. The later proof of Lin Q also uses the 
generalized energy inequality. 

To deal with localization we combine the theory of paramultiplication following the 
notes of Tao |^ . These ideas in the use of Littlewood-Paley theory were developed 
by Bony (see e.g. and Coifman and Meyer (see e.g.|Q). A good exposition can 
be found in |8|. We combine this with the theory of type (1, 1 — e) pseudodifTerential 
operators which is also described in 1^. The latter theory allows us to localize on 
cubes describing our "wavelet coefficients" . Mimicking Caffarelli-Kohn-Nirenberg's 
Proposition 2 we try to prove regularity at points not contained in any squares in 
which too much dissipation occurs. For any such point, we derive a certain level 
of regularity which we refer to as critical regularity. We find a square centered at 
such a point on which we can apply a certain barrier estimate which guarantees 
arbitrary regularity on its interior. This barrier estimate can be thought of as 
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a localized version of global existence with small data. (One can compare this 
estimate with results in ^ and in which local well posedness is established 
with hypotheses slightly more relaxed than ours.) But in order for this to work, it 
is important that there not be too many small dissipating squares on the boundary 
of our square. This combinatorial issue is an ingredient which seems not to have 
appeared before and which restricts us to the case a > 1. 



In section 2, we describe the proof of Proposition 1.1. In section 3, we describe a way 
of measuring the Hausdorff dimension of sets. In section 4, we introduce a dyadic 
model for Navier Stokes which illustrates the ideas of our proof. In section 5, we 
recall some Littlewood-Paley theory. In section 6, we describe what happens in the 
way of localization of dissipation and energy flow. Most of the general microlocal 
analysis occurs in this section. In section 7, we prove a certain minimal, what we 
call, critical regularity which is maintained away from a 5 — 4a dimensional set. In 
section 8, we describe a frequently occurring class of squares which have no highly 
dissipative square on their boundary. In section 9, we apply the results of sections 
7 and 8 to bootstrap the critical regularity into arbitrarily good regularity. 

Remark on constants and notation 

For us (, } is alway an pairing in space. 

Throughout this paper the expression A < B means A < CB where C is a constant. 
But what is a constant? Constants may depend on T, the norms of the initial 
values, and on an e which we keep fixed throughout this paper. They may not 
depend on a particular square or a particular scale which we are estimating at. 
(However sometimes we choose constants so that there is a particular initial scale 
on which constants can depend. But in that case, the point is that we make 
estimates depending on those constants for all scales). Sometimes our hypotheses 
or our definitions have constants in them. When this is the case, our conclusions 
are understood to be a multi-parameter family of conclusions depending also on 
the constants in our definitions. We will invoke propositions and lemmas using 
convenient constants in the hypotheses and we do not make any attempt to keep 
track of these. If we did, the paper might become unreadable. 
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2. Preliminaries 



We begin by recalling the proof of Proposition pTl] . 
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Proof [Proof of Proposition 1.1| 



We will use the standard notation that H'^ denotes the Sobolev space over 
with f3 derivatives. We will let W^'^ denote the Sobolev space over R.^ with (3 
derivatives. 

We first recall the energy inequality obtained by pairing the equation (|^) with u 



2 dt 



A)^«||i. <-|HI^^„ + IHIi. (3) 



From this, we obtain by integrating over time (and observing that the norm is 
always positive), that if the solution u remains smooth up to time T, we have the 
estimate 



r M%.dt<ii+T). 

Jo 



(4) 



We now pair (|l|) with {~A)>^ u in order to estimate obtain 

+ {u ■ Vu, {-A)f'u) = ~\\{~A)^u\\l,. (5) 



2 dt 

Clearly, we must estimate the nonlinear problem term 

{u ■ Vu, {-A)f^u) = ((-A)^(ii ■ Vu), (~A)^u). 
Since u is divergence free, we can bound the absolute value of the above by 

ll"llwi-p||"llvKf.<j||u|k^, 

with ^ + I = \- Now since we are assuming a > |, we obtain by the Sobolev 
imbedding theorem, 

Il"llwi-p||'"||iyf.<j|kllif5 < ll"l|ff°l|w||i/°+f*||u||H/3, 

since we must spend > | derivatives to get all three terms in and we have spent 
2a - 1. 



Applying Cauchy-Schwartz, we get immediately 

l|w||H°l|w||H= + J3||u||^fi < (5||-u||^„+^ + 

Combining this with (|^), we get 

<- II ||2 II ||2 2 

In turn, combining this with (jj) and with Gronwall's inequality gives global solv- 
ability. ■ 
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3. Hausdorff dimension 



The purpose of this section is to recall the definition of Hausdorff dimension and 
to explain a simple construction which automatically produces sets below a certain 
desired dimension. 

Given any set E C M", we define its c?-dimensional Hausdorff measure H'^{E) by 

n''(E) = lim inf V r(Bf. 

s^oceCsiE) ^ 

Here C5 is the set of all coverings of E by balls of radius less than or equal to 5, 
and if i? is a ball, r(_B) is its radius. 

We say the Hausdorff dimension of E is the infimum of the set of all d for which 
TC'^{E) is 0. Note that for the purposes of the definition of dimension it doesn't 
matter what basis of objects we cover by (balls, squares, ellipsoids of a fixed mul- 
tiplicity) as long as they are a regular family. (See [^.) 

All of our arguments are built on the notion of discretizing at a certain scale . 
We would like to be able to think of sets of a certain Hausdorff dimension as sets of 
2^'^ elements at scale 2~^. We make this intuition precise for our purposes by the 
following. 

Lemma 3.1. Let Ai, . . . , Aj, . . . be a sequence of collections of balls in R" so that 
each element of Aj has radius 2^^ . Suppose that ^{Aj) < 2^'^. Define 

E = lim sup Aj , 

j >oo 

to be the set of points in infinitely many of the UseAjB 's. Then the Hausdorff 
dimension of E is at most d. 

Proof Since the Hausdorff dimension of E is the infimum of the set of all 7 for 
which W{E) = 0, it is enough to prove that W{E) = for aU 7 > d. 

Pick j such that 2^^ < 5. We can cover E by the Ufe>j [JBeA^ B. 

Now 

■WiE) <^2'"^[2-''y, 
k>j 

which limits to zero as j goes to 00 whenever 7 > d. ■ 



4. Dyadic heuristic 



Before going into the details of the proof of Theorem 1.2, we present a heuristic to 
illustrate the main ideas. 

We shall replace our vector valued function u by a scalar valued one. We shall 
represent it by a wavelet expansion. This shall be left rather abstract. We do not 
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care which system of wavelets we use precisely so long as the size of the coefficients 
of a function with respect to the basis accurately reflect the regularity properties 
of the function. We denote our orthonormal family of wavelets {wq}, with wq the 
wavelet associated to the spatial dyadic cube Q of sidelength 2~-'('5). Then u can 
be represented as: 

Q 

where uq is the wavelet coefficient associated to the cube Q. 

Until the end of this section we reserve notation wq to denote the wavelet localized 
on a cube Q, while uq, and vq are reserved for wavelet coefficients associated to a 
cube Q. 

We need a few definitions in order to make our model precise. A cube Q in M.^ is 
called a dyadic cube if its sidelength is integer power of 2, 2', and the corners of 
the cube are on the lattice 2'Z'^. A dyadic child of a dyadic cube Q with sidelength 
2' is any of 2'^ dyadic cubes Q' contained in cube Q and such that their sidelengths 
are 2'"!. 

We will call {Q, Q' , Q") a cascade if Q' is a dyadic child of Q, and Q" is a dyadic 
child of Q'. 

Now we are ready to define an operator which will mimick the behavior of nonlinear 
term u ■ Vu. Note that we have 

IKIk~~2^. (6) 

On the other hand 

||Vi«q|U. ~2^ (7) 
since our wavelets wq 's are orthonormal and localized to the frequencies around 2^ . 

Having in mind and (0) we define a cascade down operator as follows: 

Cd{u,v)^ ^ 2'^2^'-^'>UqVqiWq". 
cascades 

Similarly we define a cascade up operator as: 

Cy,{u,v)= ^ 2'^^2^'-'^'^uq,,vq'Wq. 

cascades 

(Note that the choice of definition of the cascade down operator is somewhat arbi- 
trary. We are happy to take any bilinear operator with few coefficients which shifts 
energy to high frequency only a few scales at a time. Once cascade down is chosen, 
cascade up must be one of its adjoints.) 

Also let us define 

c(u, v) = Cd{u, v) - c„(w, v). 
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Obviously, 

{cd{u,u),u) = {cu{u,u),u), 

which impHes 

{c{u,u),u) = 0. (8) 

Having defined operator c(u, v) we can speak about dyadic version of Euler as well 
as Navier-Stokes equations. More precisely by dyadic Euler equation we mean: 

— + c{u,u) = 0. 

For the remainder of this section, we define Laplacian as A{wq) — 2^^wq. Now we 
introduce dyadic version of Navier-Stokes equation as: 

— + c{u,u) + Alt = 0. 

Also we could speak about dyadic Navier-Stokes equation with hyper-dissipation 
by which we mean: 

nil 

— + c{u,u) + {A)"u = 0. (9) 
dt 

A simple consequence of (||) is conservation of energy for all three equations, which 
is an important feature of these equations preserved in the dyadic model described 
here. 

Now we consider dyadic model for the Navier-Stokes equation with hyper-dissipation 
(H) . We would like to estimate Hausdorff dimension of the set of singular points at 
the first time of blow up, T. The nonlinear term c(u, u) on scale j looks like 2^Uq 
(if we imagine for a moment that all neighboring cubes have coefficients of roughly 
the same size), while dissipation term gives decay like 2^°'^uq. This means that as 
long as uq < 2~5(5-4a)^ ^jjg growth of uq is under control. But let us check what 
happens if uq > 2~5(5-4a)^ (This idea is made more precise by Lemma 9J. The 



reader is encouraged to work out its simpler dyadic analogue). We call this bound 
on the coefficients critical regularity and this is what we shall prove away from a 
5 — 4a dimensional set. 



We can rewrite equation ^ in terms of wavelets coefficients as follows 
dt 



{Q',Q") 



where C(^q^q',q") = -1 if (Q, Q', Q") is a cascade, where C(q^q',q") = 32 if (Q', Q", Q) 
is a cascade and vanishes otherwise. 

Having assumed uq > 2"5(^~^") for some time t and assuming that at the initial 
time i = it is much smaller, by the smoothness assumption on the initial condition, 
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we integrate ( p^ ) in time on the interval [0, T] and obtain for one of the choices of 
{Q'l Q") giving a non-vanishing coefficient: 

2^ r \uQ,UQ„\dt>2-i^'-^"\ 
Jo 

which by Cauchy-Schwartz imphes: 



^0 



The last expression can be rewritten as 



^0 



This can happen if cither 



2^^"(/ WQ'rfi)' ^ 2-^(5-^^"), (11) 







or 

2^^"(/ dt) 5 > 2-^(5-4"). (12) 
Jo 

However having in mind conservation of energy we have 

^ Q at scale j 

Thus we conclude that dO) or ( p^ ) could happen in at most < 2-'*^^^4a) cubes Q. 
Now we invoke Lemma |3.l| and conclude that the Hausdorff dimension of the set of 
points of the equation (M) at which critical regularity fails is at most 5 — Aa. 

In the remainder of the paper, we shall prove the same result for the actual Navier- 
Stokes with hyper-dissipation and we shall prove that this critical regularity implies 
stronger regularity at points away from the bad set. 

We would like to remark that the dyadic model seems quite interesting in that it 
possesses all the features of Navier Stokes which we shall use in this paper. But 
moreover the nonlinear term posesses a built-in dispersive feature - namely that 
cascade down passes energy equally to all the squares Q" which are children of Q' 
and grandchildren of Q. This dispersive feature might be useful in proving global 
solvability for the model. We believe that working out such an argument would be 
a good first step towards the Clay problem. 



5. LiTTLEWOOD PALEY THEORY AND PSEUDODIFFERENTIAL OPERATORS 



We shall use a standard Littlewood Paley partition of frequency space. That is, 
we shall have Fourier multipliers Pj (on L^(IR^)) so that their symbols Pj{^) are 
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smooth and supported in |2-' < |^| < 3(2-'), so that Pj{^) = po{2 ^x) and so that 

We define Pj = J2k=-2 ^j+k- Notice that PjPj = Pj, since Pj is the sum of all 
Littlewood-Paley projections the support of whose symbols intersects the support 
of pj(^). We can analogously define the symbols pj = X^fe=_2Pj+fc- 

The Littlewood-Paley operators satisfy the inequality 
Proposition 5.1. 

Proof Let (pj = {pj)- Then (since we are working in M^) we have 
Now Pjf = <t)j * f so by Young's inequality, we have 

which is finite since is a Schwartz function. But by the relation between and 
(^0) we have 

||</.,-|U. =2¥||0o||^., 
which proves the proposition. ■ 

The idea of Littlewood Paley theory is that Pjf is like a combination of wavelets 
supported on cubes of length and that the previous proposition is sharp only 
when a large proportion of the energy of Pjf is concentrated in a single one of 
these cubes. Hence one is led to try to localize Pjf in space on such cubes and 
this is marginally consistent with the Hcisenbcrg uncertainty principle. However we 
have some room for error in our dimension estimate since it is a closed condition. 
Therefore we fix an e > and never try to localize better than within 2~^'^^~'^\ 

Our localization is nearly perfect if we neglect negligible quantities. For the remain- 
der of the paper, whenever we are at scale j, we neglect quantities of size < 2"^°°-' 
since they will not affect our estimates. Similarly we neglect operators whose norms 
are smaller than 2"^°"-' provided they will only be applied to functions whose norms 
are < 1. The choice of 100 is arbitrary. It is large enough so that it does not affect 
the L°° norm of u. But in fact our techniques for showing quantities arc negligible 
rely on Schwartz function properties and could give arbitrary exponent with loss in 
the constant. The current approach works well principally because of conservation 
of energy, since the norm of u is certainly < 1. 

For any cube Q with sidelength greater than 2~-'(^~^^\ we define a bump function 

4>Q.j which is positive, is bounded above by 1, equals 1 on Q and is outside of 
(1 -|- 2~^'^)Q. Further, we require for each multiindex a that there is a constant Ca 
independent of Q so that 

|£'"<^Q,il <Ca2l"l^'('-^^. (13) 
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We say that any bump function which satisfies the estimates ( |13D is of type j. 

The map 4>Q,jPj acts much like a projection, and we shaU treat \\4>Q,jPj IWl^ if 
it were a wavelet coefficient. When we deal with a cube Q of sidelength 2~^^-^^^'^\ 
we shall define j{Q) — j and we shall denote \\4>Qj{Q)Pjf\\L2 by /q. Further if 
j{Q) = j, we say that Q is at level j. 

Proposition 5.2. Given f with ||/||l2 < 1, and <f> a hump function of type j we 
have that the quantity 

\\<PP,f -P,<j>P,.f\\L^ 

is negligible. 

Proof Let us define = 0i + 02, where 02(0 = X{|4|>^2j}'^(0- We have 

0P,/ - P,0P,/ = 0iP,/ - P,0lP,/ + hPjf - ^,02P,/. 

By our estimates on the derivatives of 0, we can get 

1 102(0 1 k~ = negligible, 
while, because of the Fourier transform supports, 

P,0iP,/ = 0iP,/. 

Thus the proposition is proved. ■ 



Combining propositions 5.1 and 5.2, we get the following extremely useful lemma. 
Lemma 5.3. Let (p be a bump function of type j . Then 

WpPjfllL^ < 2^ \\,pPjf\\L2 + negligible. 

Proof We neglect negligible terms. Then we have 4>Pjf = Pj(j>Pjf, by proposition 



5.2. We estimate ||Pj0P)/||L~ by proposition 5.1 



Proposition 5.4. For any cube Q, we have that for any f with \\f\\Hf ^ 1 with 
/3 < 10 that 



is negligible. 
Proof 



(1-0 -.j(Q) ^ ^ )Pj't'Q,j(Q), 
(1+2 2 )Q,j{Q)' ^'■"^^> 



is a composition of type (1, 1 — e) pseudodifferential operators whose symbols have 
disjoint support. Thus it is smoothing. ■ 

Notice that the above proposition really says that we can move bump functions 
across Littlewood Paley projections as long as the bump functions proliferate and 
increase in support. In other words, the proposition can be rewritten as 

'^(i+2'^2'-")Qj(Q)^-'"'^Q-j(Q) " Pj(t)Q,j(Q) + negligible. 
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Similarly, this can be used to remove inconvenient bump functions, provided there 
is a smaller bump function already present in the expression. 

Proposition 5.5. Let C be a covering of a set E by cubes of sidelength 2~^^'^~'^'^^ 
for some fixed j . Then for any f £ , we have 

\\XEP,f\\h<J2fQ- 
Qec 



Proof We simply observe that 



Qec Qec 



Essentially what we have done in this section is to use approximations to projections 
which are uniformly pseudodifferential operators of type (1,1 — e). The negligible 
interaction of distant squares can be just as well derived from the asymptotic for- 
mula for composition of such operators. The composition of two (1, 1 — e) operators 
whose symbols have disjoint support is infinitely smoothing and hence negligible. 



6. Dissipation, Trichotomy, and Forbidden squares 



In this section, with u a solution to the Navier Stokes with hyperdissipation, we 
write down estimates involving the growth of the mq's. This is the only section in 
which we work directly with the equation. By the end of this section, we will have 
reduced our problem to a combinatorial one about ordinary differential equations 
on graphs. This section is presented with apologies to Terry Tao, since much of it 
is based on a course he gave at UCLA in Winter 2001. 

Recall we have 

— +u-Vu + Vp=-(-A)> (14) 

together with 

V-u = (15) 
In light of (|l5|), we can rewrite (|^ as 



— +T(wVu) = -(-A)V (16) 

where T is the projection into divergence free vector fields. The operator T is a 
singular integral operator and is also a Fourier multiplier. We pick a cube Q of side 
length 2"-'*^^"^'^) and compute the pairing of the equation with PjcjjQ jPjU. We 
obtain the energy estimate 
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We shall estimate the two terms on the right hand side of ( p7| ) separately. 
As before we use the notation 

2 

fe=-2 



This has the advantage that 



P — P P 



We make some definitions. We define for each cube Q, the set Af^{Q), the nuclear 
family of Q to be a union of sets Aq,Bq,Cq,Dq, Eq, where Aq,Bq, Cq,Dq,Eq 
are covers of Q = (1 + 2^)Q by fewer than 1024 cubes each at levels respectively 
j — 2,j — + 1, and j + 2. We define recursively M^{Q) to be the union of all 
N^{Q') for aU Q' e M'-^'^{Q). Thus in particular, we have 

#{^^\Q)) < 2^^', 

since 5{2^°) < 2". 

Proposition 6.1. Let Q be a cube and j = j(Q)- Then 

((-A)"m, Pj(t)l jPju) > c22"Ju| - c'2(2"-2^)J ^ u^, - negligible 

Q'&M^iQ) (18) 

Proof Note that 

{{-l^ru,P,4>l,,Pju) = {{-Arc^Q,,P,u,^Q,,P,u) + {[cl,Q,,P,,{-Ar]u,ct>Q.jPju) = X+Y. 

Note that by proposition |]^, we have 

X ^ 22"^M^ - negligible. 

To estimate Y, we observe that [4>Q,jPj, (—A )"] is of order 2a — e since (pQ.jPj is 
of type (1, 1 — e). Further by proposition |5.2| , we have 

[0Q,jP, , (-A)"]P, - [0Q,,P, , (-A)"] = negligible. 



Even further, applying proposition 5.4, we get 



^ = (['/'Qj-P)> (-^)"]'/'q,j-F'jW: 'PqjPj'^) + negligible. 

Now applying the mapping properties of operators of order 2a — e and Cauchy 
Schwartz, we get for any number K , 



\Y\ < l22"^-«^ + X2(2"-2^)^||0Q/,^||i. 



K 

Now applying proposition |5.4 we get the desired result. 
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Let US define for every I, 

2 _ \ ^ 2 



Then we have 

Corollary 6.2. Let Q be a cube and j = j{Q). Then for any I, 

J2 ii-^ru, P,4>Q'^,P3U) > c22"^uX',(Q) - C2(2"-2^)^«^,^,(^^ _ neghgible 
Q'eM'iQ) (19) 



Proof Simply sum proposition 6.1 over M (Q) 



We want to use Proposition |6.1| to fight against the growth of uq. Thus the term 



2(2a-2e)jy|^ WOuld appear to be a serious nuisance. However we are able to show 

Lemma 6.3. For any interval of time J C [0,r] and any cube Q we may find an 
so that 

J ujj-i^Q^dt + negligible>2^''^ J u^j-i+i^Q-^dt. (20) 



Proof By conservation of energy, for every I < ^ and every t we have, 

■"a/'' + i(Q) ^ 1- 

Since our constants can depend on T this means 



Now suppose the lemma is false. Applying the opposite of (^), for Z = 1 
we get 

J '"a/'2(Q)^^ + negligible < 2^^^^-' = negligible, 



400 



which implies that (20) holds for Z = 1 



We will be able to apply Lemma 6.3 together with Corollary 6.2 to show that for 
any cube and any interval in time, there is an Ith iterated nuclear family for small 
I which is undergoing dissipation. 

Now we turn our attention to the term 
We rewrite it as 

Gq = {-<t>Q^PjTP,{u ■ VU), <t>Q,,PjU). 

Now we use the "trichotomy" . We can write 

Pj[u ■ Vm) = Hj^ih + Hj^hi + Hjjih + Hioc, 
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where the low-high part is given by 
the high-low part is given by 

k<j-1220 

the high-high part is given 

HjMh^ J2 PjiiPku) ■ PkVu) + PjiiPku)- PkVu),, 

(technically that was the end of the trichotomy) and the local part is given by 
Hloc^ Pj{{Pku)-Pk^u)+ Y Pj{{Pku)-Pk^u),. 

Now we break up Gq in the obvious way: 

Gq = Gqjh + Gq^m + Gqmi + Gq^ioc, 

where 

GQ,lh = {-(f>Q,jPjTHj^lh,(f>Q.jPjU). 
Gqm = {-(j)Q,jPjTHjMi(t'Q,]P]U). 

Gq.hh = {-(f'Q.jPjTHj^hh'pQjPju)- 

Gq,Ioc = {-(f>Q,jPjTHj^loc,(t)Q,jPjU). 

Now fixing Q and j we would like to estimate each of GQ^ih,GQ^hi,GQ^hh, and 

We take a moment to be careful about how we localize. For any k < j, define 
Qfc = 2(J-'=)(i-2e)Q^ Pqj. j^jjy ^ > J (including j) define Qk = (1 + 2^)Q. We 
retain this notation only for this section and the following one. 

Then we have 

Lemma 6.4. For any 5 > 

j 1000 

\GQ,ih\ + \Gqm\ < Y 2^^^MQfcUA^i(Q)WQ + 2^^^+^^uj^i(Q)UQ + negligible. 

k=Sj 

Proof We consider for k < j — , the relevant expression 
Gqjh.k = {-(t)QjPjT{PkuPj\/u), (l)QjPju). 

We divide this into two cases which are k < 5j and k > dj. In the second case, we 
use the idea of proposition (we just have to use a higher degree of smoothing 
than is used in that proposition) to observe 

GQ,ih,k = {~(l}QjPjT{{(t)Q^,kPk)u ■ V((/)g^. jPju)), (pQ^jPju) + negligible. 



CAFFARELLI-KOHN-NIRENBERG 



15 



(This is because V acts on 



only where it has a neghgible effect on the whole 



quantity by proposition 5.4.) Now we simply observe by proposition 5.3 that 
\\4'Qk,kPku\\L^ < 2^uq^ + negligible, 



and by proposition 5.4, the proof of proposition |5.2| (to control the action of V) 
and proposition |5.5| that 



(Q) 



negligible. 



On the other hand for the first case, there is no point in localizing to Qk because 
k is too small and so the error is too big. Thus in in this case, we simply estimate 



\PkU\ 



35j 

< 2 — u, 



Qk- 



Summing these estimates gives the desired bound for Gg^ih- The bound for Gqm 
proceeds likewise (and gives a better estimate since the derivative falls on the level 
k term.) 



Lemma 6.5. 



\Gq,ioc\ < 2 2 MQM 1000 + negligible. 

JV e 



Proof From the definition of Gq,ioc we have: 

2 3 + ^ 

Gq.ioc^Y.^ J2 {<^Q,jP3T{Pk+iu-^Pku),ci>Q.,Pju)+ {'l>Q.jPoT{PkU-VPk+iu),^Q,jP,u)]. 

fc=j_1000 



Since the above sum has only < 1 many terms by applying Proposition 5.4, we 
observe that for some particular values of k,l 

\Gq.Ioc\ < \{(f>Q.jPjT{(j)Q.jPkU ■ (t)Q^^jVPk+iu),Pj(t)Q,ju)\ + negligible. 

Now we apply proposition 5_^ together with Cauchy Schwartz to obtain 



3j 

ll'/'Qj ^ 2~||0Q^ jPfcw||i2 + negligible. 



Using proposition 5.5, we observe that 

Finally direct calculation shows 

\\(t^Qj,3'^Pk+iu\\L'^ ^ \ \^<I>^^^^^^Q Pk+iu\\ < 2^ Uj^icm^^^+ negligible, 

where the first inequality comes from the fact that 6 -2ej = 1 on the support 

(1+2 — )Q.j 

of 4>Qj ,j and the second inequality comes again from proposition p.q . Combining 
all these estimates proves the lemma. 



16 



NETS HAWK KATZ AND NATASA PAVLOVIC 



Lemma 6.6. 

\Gqmh\< E "Q2¥+'=||</.Q^^,P,u||i.. 



Proof We begin similarly to before by estimating 

{(t^Q,jPjT{PkU ■ WPku),(f>Q,jPju). 

With the other terms one can proceed likewise. Now, similarly to what we have 
already done, we observe that we can write this as, 

{(I^QjPjTi'l'QjjPkU ■ W(l)Q^jPku), <j)Q^jPju) + negligible. 

Now we estimate 

\\4>q.jPju\\l'^ < 2~uq + negligible, 

2 

IIWQj j-Pfcw||L2 < S'' ^ \\(f)Q.,jPk+iu\\L2 + negligible. 
Combining these estimates gives the desired inequality. ■ 



In fact, we get a slightly better estimates, for instance by the div-curl lemma, but 
it does not seem to be necessary. 

The previous three lemmas are somewhat wasteful, particularly in the definition of 
the local part. We are using more squares than we need to cover. For the result of 
the final section we need a somewhat more efficient decomposition, which is proved 
in exactly the same way. We formulate it here. 

For Q a cube at level j, we define the ancestors of Q, 

-^(Q) = {Qk}e]<k<j-4:- 

For fc < j — 4, we define the collection Sk{Q) to be a covering (with overlap < 1) 
of the cube Qj-4 by cubes at level k. We define the strict extended family 

j I 1000 

and we define !F{Q), the followers of Q by 

HQ) = u ^'^(Qy 

k>j+1020 



Then we have the following corollary of the proofs (but not the statements) of 
lemmas 3.4, 3.5 and 6.6 

Corollary 6.7. For any 5 > we can estimate 

Gq < Zq^ihhi + Zq^oc + Zq^hh + i^^^^^i^^qu^^Q), 

where 

Zq,ihhi = ^ ^ 2^+"^uq^uqruq, 

k<3-iQ'&N'^{Q) 
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2 



and 

Zq,Ioc = 2^ ^ UQUQ'UQn, 

Q',Q"eS{Q) 

and finally 



Now we are ready to describe our singular set. We will say that a cube Q of 
sidelength 2-^(1-2^) is bad if 

r [ V22"'=|0Q,,F,u|2 > 2-(5-4").-ioo.. . (21) 

Let Ej be the union of 2^^Q for all cubes Q of sidelength 2^^^^~'^'^^ which are bad. 

We will need the following well known covering lemma of Vitali (see |^): 

Lemma 6.8. Let C be any collection of cubes, then there is a subcollection C so 
that any two cubes in C are pairwise disjoint and so that 

\JQc U 5Q. 
Qec Qec 

Proposition 6.9. There is covering Qj of Ej by cubes of sidelength 2^^^^"'^'^^ so 
that 

Proof Let C be the collection of cubes 2Q where Q is a bad cube at level j. From 



Lemma 3.5 we know that there are disjoint cubes {2Qp}p^z such that the collection 
consisting of {lOQp} covers the set Ej. Any cube of sidelength 10 • 2~^^^~^'^'> can 
be covered by 1000 cubes of sidelength 2~^^^~^^' (and 1000 is a constant.) We will 
define Qj to be the covering formed by the union of the thousandths of the elements 
of {lOQp} . 

In order to count cubes in Qj it is enough to count the disjoint cubes used in the 
construction of Qj. 

However we know since the 2Qp's are disjoint that 

r22"'=|0Q^,,P,z,|2 < r f |A"Kp < 1 (22) 
p k -^0 Jo J 

by conservation of energy, while we know 



p k 

by the badness of the cubes 



„ Jo 



(5-4Q)j-100ei 



(23) 



Combining the inequalities ( p2[ ) and ( p^ ) implies the claim. 
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We pause for a moment to apply Lemma p.l[ 

Corollary 6.10. The dimension of E — limsup^ ,00 hounded by 5 — Aa 

0{e). 



If we could show that if a; ^ i? then a: is a regular point of u at time T - in other 
words that 

limsup t)| < 00. 
t — >T 

(Indeed if we could show it is a regular point for a derivative of u of any fixed order) 



this would immediately imply Theorem 1.2 



Unfortunately, we will show something slightly more complicated - the same state- 
ment for a somewhat larger collection Q'^. 

In the following section, we will investigate what are the immediate consequences 
oi X ^ E. Saying that x ^ E is the same as saying that there exists a j so that 
so that for any k > j, we have x ^ Ek- Denote Fj as the set of points with this 
property. To prove a regularity statement about E, it suffices to show that that 
statement holds for any jo, provided x € Fj^. However fixing j, we may change 
our constants in the definition of bad square so that x is not contained in any bad 
squares. Thus we may as well assume x is contained in no bad squares. This will 
be our hypothesis in the following section. (And as promised, the constants will 
now depend on jo-) 



7. Critical regularity 



Heuristically, the worst part of Gq is Gq^ioc which looks at scale j like 2^^ Uq. On 
the other hand, dissipation gives decay like 2^"-'uq. This should mean that as long 

as Uq < 2^'"~3)-', the growth of uq ought to be under control. It is this estimate 
which we will show in this section for any Q not contained in Ek for any k > jo- 

Theorem 7.1. Let Q be as above, then there is a constant C , depending only on 
T , the initial conditions for u, the constant in the definition of badness, and jo so 
that 

UQ{t) < C22("-t-i)jW). 

Proof We proceed by contradiction. Suppose the theorem is false. We let Tq be 
the first time and Q be the largest cube so that 



Uq 



{To) > C22("-t-i) 



Now since our initial data is smooth, at the initial time, (since j is chosen sufficiently 
large), we have 

uq{0)<2-'°°°^. 
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Therefore, it must be the case that 

(■To 



I \GQ{t) - c2^"'ul + C2(2"~2^)J ^l')dt > 24(" 



By using proposition 6.3, we can replace Q by an extended nuclear family A/''((5), 
with / < ^ for which 



For the current theorem, this is all we need. Since Uj^i^^q^ also begins < 2 lo^OJ^ 
we must have 







(Q)' 



where we define 

GMi(Q){t)^ ^ Gq^t). 

Since there are only < 1 cubes in N^{Q)^ for this to be the case, there must be a 
Q so that 



To _ i-T 



We contradict this by using Lemmas |6^, 6.5, and 3.6 to estimate Gq. By our 
definition of Ej^ the cube Q is contained in no larger bad squares. 

Each estimate contains a factor of Uq which we take out using the estimate \uq\ < 
22(a-|-3)J j which we get from the definition of Tq. 



First let Qi be a nuclear family member of Q and Q2 be a distant ancestor at level 



k. Suppose k < dj, we see that 2*^^+^2^)-'mqj^ < Thus we need not worry 



about this term. 

Suppose k > 6j. We must estimate 

fT 



2 '2 +^UQ,UQ^dt. 







Note that both Qi and Q2 are good squares. Thus we have the estimates 







and 

< 2(2"-3)j-10ei 

Qi — 
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Applying Cauchy-Schwartz, we get (using a > 1) 
Jo 

Summing over (there are only j terms) provides the desired estimate on Gq + 
Gq if^. We can estimate Gq in the same way (by allowing k as large as j + i^^.) 

We are left to estimate Gq We fix a scale k > j and pick fc2 within 2 of k. Now 
we are left to estimate 

/•To 

/ 2^ + '^\\4>Q,Pku\\L2\\^Q,^Pk,u\\L-. 

Jo 

However since Q is a good square, we have 
Thus 

J 2'=+¥|0Q ,^F,up < 2(^W-i"^J--(2"-i)('=-J), 

which is an estimate that decays geometrically in k when a > ^. By using the 
similar estimate for applying Cauchy Schwartz and summing over fc, we get the 
desired result. ■ 



8. Combinatorics 

In this section, we begin with the set Ej which is the union of a collection Qj of 
cubes with sidelength 2~^^~'^'^^^ having cardinality < 2(^~'*"+^°'''^)-'. (We assume 
Aa-A> 200e.) 

Theorem 8.1. There exists a sequence of collections Q'j of cubes of sidelength 
2-(i-2e)j #iQj) < 2(5-4"+i00£)j^ so that for any Q of length 2-(i-2e)j ^f^^^f^ 

does not intersect any element Q'j there exists a number ^ < r < 1 with the following 
property: For no k > j is there Q G Qk so that 

lOOQ n dirQ) / (24) 

Proof We refer to the elements of Qj as the bad cubes. We say that a cube of 
sidelength 2~^^~'^'^^^ is very bad if either it intersects a bad cube of the same length 
or it intersects more than ^2*^^"^"+-'^^°'^'^'^"^^ elements of lOOQfc for some k > j with 
c a small constant to be specified later. Let E'j be the union of all very bad cubes 
of length 2~'^^~2'^^ , then by the estimates on the cardinality of the Qfe's and by 
the Vitali lemma, we can see that that E'j can be covered by < 2(5-4"+ioo^)J cubes 
of length 2^*^^^2'^)-'. We refer to these cubes as Q'j,. Now we need only prove (p^. 

Let Q be a cube of length 2"'^^^'^'^^^ which does not intersect E'j. Let 'D''{Q) be the 
set of elements of lOOQk which intersect Q. Then we have the estimate 

#{V''{Q)) < c2^5-4a+150eKk-j) ^ 
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Let /fe(r) be the function defined from 5 to 1 which counts how many elements of 
T>'^{Q) intersect d{rQ). For each Q' S 'D'^iQ) define rg/ to be that number so that 
the center of Q' lies on d{rQiQ). Then 

fk{r) < ^ X(rQ,-100(2(i-2«)O-fe)),rg, + 100(2(i-2«)O-fc)))- 

Thus 

ll/fcWIki < 200c2(4"-4-i52OO-fc)^ 

Since 4a — 4 > 200e, this estimate decays geometrically with k. By choosing c 
sufficiently small, we may arrange that 

k>j 

Thus by Tchcbychcv's inequality, there must bo a value of r between ^ and 1 so 
that fk{r) = for all k > j. This is the value of r that we choose. ■ 



9. Barrier estimate 



In this section, we prove regularity on the interior of a cube Q, provided that one 
has critical regularity for cubes containing it and cubes Q' for which dQnlOOQ' 7^ 0. 

If Q is a cube and Qi C Q, we define d{Qi), the graph distance of Qi to the 
boundary of Q by d{Qi) = fc — 1, where k is the smallest positive integer so that 

2^Q^ ndQ^fH. 

Lemma 9.1. Let Q be a cube. Suppose we know that for all t < T we have that 
for any cube Q' so that Q C Q' with sidelength of Q' being 2~'(i~2'^), we have that 

\uQ,{t)\<2~'('^\ 

and suppose further that for any Q' so that lOOQ' n dQ ^ with sidelength of Q' 
being 2~'(^~2^) and I > j — 2, we have that 

\uQ'it)\<2-'('-=^\ 

then for any Qi C Q of length 2~'=(^~^*), we have the estimate 

I^Qx(i)l<2-'=^(«^\ (25) 

where 

= min(10, + — ). 

Proof At time t = Q for all Qi we have 

< o-lOOOfe 

when Qi has sidelength 2~(i~2'^)'=. Thus at time 0, the lemma is satisfied. Let to 
be the first time at which the lemma fails and Q\ be one of the cubes for which it 
fails. It must be the case by hypothesis that 32(5i does not intersect dQ. Then we 
have 

UQ,(io)~2-MQi). 
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Let ti < to, be the last time before to when 



Then we have 

j^{ulJ>2-''^^^^\ (26) 
However on the time interval (ti,to) the lemma is satisfied. 



We will invoke Corollary 6.7. Now for any Q2 G £{Qi) we have 

diQ2) > d{Qi) - 5. (27) 

This is because Q2 C 32Qi. Further for any ancestor Q3 £ A{Qi) with Q3 having 
sidelength 2-(i-2e)'^ ^ith eA: < / < - 4, 

d{Q3)>d{Qi)~5{k-l). (28) 

Further for any follower Q4 e J'^{Qi) (that is a cube which contributes to Gq^^hh 
and is in particular contained in |(5i) of Qi with (54 having sidelength 2^^^^^'^)' 
with I > k + we have 

d(Q4)>d(Qi) + (^). (29) 
Applying (p9|), we conclude that p{Qi) — 10. Thus we conclude 



ZQ,,hh < J2 2^23('-'=)2-20/2(-P(Qi)+T%)fc^ 

l>k 

by counting the elements of Si{Q) by 2'^('~''). Calculating, we find 
which since p{Qi) < 10, cannot possibly account for (E9). 



Now applying ( p7[ ) to Proposition 3.1, we observe that for any t between ti and 
to, we have dissipation at Qi of > 22"'=M^^(t) > 2(2"-2p(Qi)+f Thus to reach a 
contradiction, it suffices to show that on this time interval 

For this reason we can also ignore the ultra-low term 2-''-^+^*-'wqj^u_/v'i(Qi) in Corol- 
lary EtI 



Using (p7|), we observe that for any Q2 in £{Qi) we have the estimates 
as well as 

(4Q-5~2e)fc 

< 2 

by the lower bound on p. Now we calculate 

Q2,Q2e£(Qi) 

Thus Zq-^^ioc cannot contribute to the growth. 
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Now to estimate Zq^jiiih, we observe that for Q[ e M^{Qi) we have the estimate 

un' < 2("''^'5i)+to)'= 



while for ancestor Q4 of sidelength 2 ^"^^'j we apply ( p8| ) (as well as the hypothe- 
ses of the lemma for squares larger than Q) 



< 2 ''*°"2"''"' +(j-100e)(fc-i) 



since of course | > ^ „'^° + 200e. Now we just estimate 



Thus could not have grown which is a contradiction. 



But now we we have in fact proven the main theorem. We define Ej as in section 
|[ We need to show that u is regular at any x not contained in any E'j with j 
larger than some integer j^- By changing our constants, we can say x is not in in 



any E',. Now by theorems 7.1 and B.l, for any cube Q centered at x, we can find 



a cube almost half as large which satisfies the hypotheses of Lemma But the 
conclusion of the lemma implies that x is a regular point. Thus any singular point 
must be contained in E = limsup£^'. By lemma 3.1, we have dim(i?) < 5— 4a+20e. 



Letting e tend towards 0, we get Theorem 1.2. 
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